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Ein verbessertes Schätzverfahren für Gesamtpopulationsgrößen 
bei Tagfaltern und anderen Invertebraten 


M.A. PFEIFER 


Abstract: Anadvanced method for estimating total population sizes of butterflies 
and other invertebrates. — During the phenology of most species of butterflies all 
individuals of a population are never present at the same time. For this reason they are 
different from many other species frequently analysed by mark-recapture methods. 
Therefore a method for estimating butterfly populations sizes has to be fitted to this 
specific demand and an adequate method for estimating total population sizes based on 
the Jolly-Seber method has been developed. 


With the Jolly-Seber method the population size N; at time of sample 1 is estimated. 
Additional estimations are given for the number of new animals B; joining the 
population in the interval i to 1+] and which are still alive at time of sample i+1. Also 
an estimation for the residence time 6, between i and i+} is delivered by the Jolly- 


Seber method. Nevertheless, the Jolly-Seber method does not meet the problem of 
estimating population sizes of species in which not all individuals belonging to a 
population are present at one point of time, like most butterflies do. 


SETTELE et al. (2000) introduced a procedure based on the Jolly-Seber method for 

estimating the total population size adding ali joining individuals over the entire time of 
k-2 

the sampling period: Nj, = N2 + > B; . But this procedure does not take account for 
i=2 

the number of joining individuals between two samples i and 1+] which already have 

left the population at the time of sample i+1 (either they died or emigrated between 

samples 1 and i+1). 


Therefore a correction factor for B; has been established. In the case that individuals 
joining the population at definite times — e. g. in butterflies hatching at a definite time 


sng Ni l-o? 

of the day - Bi has to be adjusted by multiplying with Kl; I 
number of days between sample i and i+1; t is the span of time (as a fractional amount 
of a day) between the definite time when individuals enter the population (e. g. time of 
hatching) and time of sample i+1. 


-lnġ; . : ; : ai 
K2 sang, is the factor of correction for species, which enter the population 
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continuous, or for the case that the time ıs unknown when individuals enier the 
population. 


To minimize the effort for estimating total population size one would like to extent the 
span of time between samples, but with extending the time span the precision of the 
estimation will be reduced. Therefore the span of time between samples ıs always a 
compromise between precision and expenditure of human work. Nevertheless some 
rules of tomb should be considered: In small populations and in populations with a low 
daily residence time the time span between samples should be short. 


Key words: Population size estimation, total population size, Jolly-Seber, 
invertebrates, insects, butterflies. 


Einleitung 


Populationsgrößen werden weitgehend als momentane Anzahl von Individuen verstan- 
den. Bei Tagfaltern kann die Festlegung auf die momentane Populationsgröße meistens 
keine Anwendung finden, da zu keinem Zeitpunkt alle Individuen des fliegenden Stadi- 
ums (= Falterstadium) während einer Flugperiode gleichzeitig adult sind. Exemplarisch 
sei die Bläulingsart Maculinea nausithous genannt. Sie ist univoltin mit einer Flugzeit 
regional verschieden beginnend zwischen Mitte Juni und Mitte Juli, welche sich über 
mehrere Wochen erstreckt. Die mittlere Verweildauer beträgt in den bisher untersuchten 
Populationen nur | bis 5 Tage (PFEIFER et al. 2000). Zu keinem Zeitpunkt sind alle Indi- 
viduen einer Population gleichzeitig adult. Somit kann mit der Jolly-Seber-Methode 
nicht direkt die Anzahl der Tagfalterimagines geschätzt werden, die über die gesamte 
Flugzeit hinweg der Population angehörten. 


Bei Tagfaltern interessiert also vielmehr die Anzahl aller Imagines innerhalb einer Flug- 
periode (SETTELE et al. 2000, MUNGUIRA & MARTIN 1999). Populationsgrößenschätz- 
methoden sind bisher nicht auf dieses spezifische Problem ausgerichtet. Sie schätzen 
lediglich momentane Populationsgrößen. 


SETTELE et al. (1998; 2000, Corrigenta unter http://www. .ufz.de/spb/nat/ settele/ tagfalter. 
html) nennen eine Methode zur Schätzung der so genannten Gesamtpopulationsgröße. 
Von BIEDERMANN (1997) wurde dieselbe Methode zur Schätzung von Gesamtpopu- 
lationsgrößen an Zikaden angewendet. Diese Methode führt allerdings in nahezu allen 
Fällen zur Unterschätzung der Populationsgröße, wie nachfolgend erklärt wird. Es wer- 
den Korrekturverfahren vorgeschlagen, welche diesen Bias weitgehend kompensieren. 


Bisherige Methode zur Gesamtpopulationsgrößenschätzung 


Tagfalterpopulationen sind offene Populationen. Als offene Populationen werden Popu- 
lationen bezeichnet, welche während des Untersuchungszeitraums Zu- und Abgänge 
aufweisen. Zugänge kommen bei Tagfalterpopulationen durch Einwanderung und frisch 
geschlüpfte Individuen zustatten. Abgänge aus der Population beruhen auf Emigration 
und Tod von Individuen. 


Größen offener Populationen lassen sich mit Hilfe der Jolly-Seber-Methode (z. B. SEBER 
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1982, KREBS 1999, SETTELE et al. 2000) schätzen (siehe Tab. 1 für die Erklärung der 
verwendeten Parameter). Die Methode geht von der Annahme aus, dass die Verweil- 
wahrscheinlichkeiten aller gleichzeitig in der Population anwesender Individuen gleich 
groß und somit unabhängig von individuellen Faktoren, z. B. dem Alter, sind. Bei 
Tagfalterarten kann davon ausgegangen werden, dass die Überlebenswahrscheinlich- 
keiten aller gleichzeitig anwesenden Individuen unabhängig vom Alter sind (PFEIFER et 
al. 2000). Ebenso muss gelten, dass es keine altersabhängigen Emigrationsraten gibt. 
Wenn davon nicht ausgegangen werden kann, sollte zumindest eine niedrige Emigra- 
tionsrate vorliegen. Eventuell gibt es geschlechtsspezifische Unterschiede in der Ver- 
weilwahrscheinlichkeit und in diesem Fall müssen die Geschlechter selbstverständlich 
getrennt betrachtet werden. 


Tab. 1: Bedeutung der Variablen. Bei einem Verlust von Individuen zwischen zwei Fangperioden 
kann nicht unterschieden werden, ob der Verlust durch Tod oder Abwandem entsteht. Dennoch 
wird in der Literatur meistens der nicht völlig korrekte Begriff Überlebenswahrscheinlichkeit (=1- 
Sterbewahrscheinlichkeit) benutzt. Korrekter wäre Verweilwahrscheinlichkeit. Zum Verständnis 
der Variablen trägt zusätzlich Abb. 1 bei. 

Unglücklicherweise wurde als Index von B (der zum Zeitpunkt 1+1 geschätzten Anzahl neu zur 
Population hinzugekommenen Individuen) 1 in B; etabliert. Eindeutiger wäre B;+ gewesen. 


Zahl der Individuen in der Population zum Erfassungszeitpunkt i. 
Gesamtpopulationsgröße. 
Verweil-(Uberlebens-)wahrscheinlichkeit zwischen i und (i+1)-ter Erfassung. 


Tagesverweil-(Uberlebens-)wahrscheinlichkeit. 


Zwischen i und (i+1)-ter Erfassung neu hinzugekommene Individuen, die zum Zeitpunkt i+] 
noch anwesend (lebend) sind. 


Zwischen i und (itl)-ter Erfassung neu hinzugekommene Individuen inkl. der bis zum 
Zeitpunkt i+1 wieder ausgeschiedenen Individuen 


Anzahl Tage vom Erfassungszeitpunkt i bis zur nächsten Erfassung zum Zeitpunkt i+1. 


Zeitspanne (Teil eines Tages) zwischen letzten Eintritt (Schlupf) von Individuen vor dem 
Erfassungszeitpunkt von i+1 und dem Erfassungszeitpunkt i+1 selbst. 

Anzahl zum Eintrittsereignis (Schlupfzeitpunkt) j zwischen den Erfassungsereignissen 1 und 
it] neu hinzukommende Tiere. 





Verweil-(Uberlebens-)wahrscheinlichkeit zwischen den Eintrittsereignissen j und j+1. 


Die Populationsgröße N zum Erfassungszeitpunkt 1+1 errechnet sich 
N =N.® +B; (1) 


bzw. errechnet sich die Anzah] zum Erfassungszeitpunkt 1+] neu hinzugekommener 
Individuen nach 


BaNui-Ni*9; (2). 
Die Methode zur Schätzung von Gesamtpopulationsgrößen nach SETTELE et al. (2000) 
beruht nun darauf, alle neu hinzugekommenen Individuen B; aufzuaddieren: 


k-2 
Na=N +ò B © 


i=2 


116 


Die Populationsgröße bei der erste Erfassung und der Zuwachs bei der letzten Erfassung 
lassen sıch mit der Jolly-Seber-Methode nicht berechnen, weswegen diese Randdaten 
wegfallen und es daher zu einer geringfügigen Unterschätzung der Populationsgröße 
kommen kann. 


Vorrausetzung, um mit der bisherigen Methode nach SETTELE et al. (2000) die Anzahl 
aller Individuen in einer Flugperiode zu schätzen, ist, dass die zwischen zwei 
Erfassungszeitpunkten 1 und ı+l neu hinzukommenden Individuen zum 
Erfassungszeitpunkt it] noch alle anwesend sind, wie nachfolgend erklärt wird. 


Kritik an der bisherigen Methode 


Angenommen, alle in die Population eintretenden frisch geschlüpften Individuen einer 
Schmetterlingsart erscheinen zum selben Zeitpunkt. Einwanderung von Individuen spiele 
nur eine untergeordnete Rolle. Die Aufnahme der Stichprobe für die Schätzung der 
Tagfalterpopulationsgröße finde direkt im Anschluss an das Eintretungsereignis statt. In 
diesem Falle sei die Annahme erlaubt, dass zum Erfassungszeitpunkt weitgehend alle 
frisch geschlüpften und nun in die Population hinzukommenden Individuen anwesend 
sind und somit in die Stichprobe eingehen können. 


Wird somit täglich erfasst, dürften bei dem gerade genannten Beispiel nahezu alle Indi- 
viduen die Möglichkeit gehabt haben, in die Stichprobe einzugehen. 


Wird allerdings nicht täglich erfasst, wird ein Teil der neu ın die Population eintretenden 
Individuen nicht bis zum nächsten Erfassungszeitpunkt überleben und in der Population 
verbleiben. Ein Teil der zwischen zwei Erfassungszeitpunkten neu hinzukommenden 
Individuen wird also nie zu einem Erfassungszeitpunkt anwesend sein und kann somit 
nicht in die Stichprobe eingehen. 


In eine Population von Tagfaltern können während der Flugzeit täglich neue Individuen 
eintreten. Zunächst soll zum leichteren Verständnis der Fall betrachtet werden, dass 
Individuen zu definierten Zeitpunkten ın die Population eintreten und dazwischen 
Intervalle ohne Zuwachs liegen. Praktisch gesagt, also täglich Individuen zu einer festen 
Uhrzeit in die Population eintreten und bis zum nächsten Tag keine weiteren Individuen 
mehr hinzukommen. 


Die Zahl B; der zum Erfassungszeitpunkt i+1 neu hinzukommenden noch in der Popula- 
tion anwesenden Tiere errechnet sich dann folgendermaßen: 


B=((...((Bz1*z1+Bz2)* P 22+ Bz3)* Pr... +Bz,1)*Pz,1tBz,)*®z, (4) 
Wobei Z, bis Z,., die Tage zwischen zwei Erfassungszeitpunkten sind. Z, ist schließlich 
der Tag, an dem die (i+1)-te Erfassung stattfindet. 


Damit das Berechnungsverfahren zur Gesamtpopulationsgröße nach SETTELE et al. 
(2000) korrekt ist, müssen alle Individuen, die in der Zeit zwischen i und i+] neu in die 
Population eintreten, zum Erfassungszeitpunkt 1+1 noch anwesend sein, um in die Stich- 
probe eingehen zu können. Das ist der Fall, wenn keine Mortalität vorliegt 


y 
(®21=®2=...=®,,1=®2,=1). Dann gilt: B; = $ B; . 


Fl 
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Auch ist das Berechnungsverfahren bei einer Überlebenswahrscheinlichkeit ®,,<1 rich- 
tig, wenn nur zum Zeitpunkt des letzten Eintrittsereignis (entspricht hier dem Erfas- 
sungszeitpunkt von itl) Individuen in die Population eintreten (somit Bz,=Bz.=... 
=B;,.=0 und nur Bz, > 0). Zudem muss dieser Zeitpunkt des Eintretens in die Population 
gleich oder zumindest so zeitnah dem Erfassungszeitpunkt sein, dass Mortalität zwischen 
Eintreten in die Population und Erfassung vernachlässigt werden kann. 


Diese Bedingungen dürften in den seltensten Fällen erfüllt werden. In der Regel werden 
an allen Tagen Individuen neu ın die Population eintreten und an allen Tagen die Ver- 
weilwahrscheinlichkeit ®,; < 1 sein. Dadurch wird ein Teil der zwischen i und i+] neu 
hinzukommenden Individuen nicht beim Stichprobenahmeereignis i+] anwesend sein. 
Dieser Anteil an Individuen wird durch das Schätzverfahren für Gesamtpopulations- 
größen nach SETTELE et al. (2000) unterschätzt. 


Verbesserung des Berechnungsverfahrens 


Vereinfachende Annahmen 

Es gilt für den Zuwachs B; zum Erfassungszeitpunkt i+1 folgende Gleichung: 
B=((...((Bzi*®zıtBa)*®ztBz3)®P23t...+Bz,.1)*P2,.1tBz,)*®z, (4) 

Aber Bz, und #7; für je {1;2;...y} sind nicht bekannt und Aussagen über Bz; bzw. #7; 
können nicht direkt getroffen werden. Die vereinfachende Annahme sei daher erlaubt, 
dass an allen Tagen zwischen i und i+] die gleiche Anzahl an Falter schlüpfen 
(Bzı=Bz;=...=Bz,) und die Mortalitäts- und Emigrationsraten (Verweilwahrscheinlich- 
keiten) über alle Tage zwischen i und it] gleich bleiben (#,,=® ,,=...=®,,). 


Nachfolgend werden die Zuwachsraten Bz; und Verweilwahrscheinlichkeiten ®;, für alle 
je{1;2;..;.y} kurz Bz und #, genannt, da sie nach der getroffenen Annahme für jedes 
Eintrittsereignis j den gleichen Betrag haben. 


Dann lässt sich Gleichung (4) auch folgendermaßen schreiben: 
B=((...((Bz*®z+Bz)*Pz+Bz)Pz+...+Bz)*Pz)tBz)*P, (5) 


Gleichung unter idealen Erfassungsbedingungen 


Unter idealen Vorraussetzungen folgt am Tag der Erfassung Z, der Erfassungszeitpunkt 
von i+] unmittelbar auf das Eintreten neuer Individuen in die Population. Unter diesen 
Vorraussetzungen kann annähernd angenommen werden, dass am Tag y zwischen dem 
Eintreten neuer Individuen in die Population und der Erfassung keine Mortalität oder 
Emigration stattfindet. Die Verweilwahrscheinlichkeit beträgt also zwischen dem 
Eintreten der Individuen in die Population am letzten Tag und der Erfassung 1. 
Gleichung 5 ist daher für diesen Fall zu verändern: 


B=((...((Bz*PztBz)*P®ztBz)®z+...+Bz)*Pz)+Bz (6) 
Gleichung (6) wird aufgelöst: 
B, = B,*6%'+B,*6%? +...+ B, *6, +B,*6, (7) 


und vereinfacht dargestellt: 
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B,=B,*>'9) (8) 


j=0 
Für $7<1 liegt eine Summenfolge mit folgender Eigenschaft vor (geometrische Reihe): 


I-b, 


Z 





B, =B,* (9) 


Da =7,*z,_\*..*82, und angenommen wird $7;=7,=...=®z,, gilt 
>; T % (10) 


bzw. aufgelöst nach &;: 


1 
6, =67 (11) 
Bz ist die gesuchte Größe der täglich neu hinzukommenden Individuen und errechnet 
sich nach Umstellen von Gleichung (9): 


zZ 


B, = B; tz (12) 


Durch Einsetzen von Gleichung (11) ergibt sich: 


I 
1-6’ 
B,=B,*-—* (13) 
1-6; 
Somit ergibt sich als Anzahl insgesamt zur Population zwischen i und i+ł - also in einer 
Zeitspanne von y Tagen - neu hinzukommende Individuen: 


Bi Ges=Bz*y (14) 
Die Gleichung nach SETTELE et al. (1998, 2000) ist daher folgendermaßen zu erweitern: 


— AY 


] : 
, = B,*——+* 15 
iGes i 1-6, y ( ) 


1-6? 


Der Term *y wird nachfolgend Korrekturfaktor Ku; genannt. Somit gilt 





Bices = B;*Ku.. (16). 


Gleichung unter pessimistischen Erfassungsbedingungen 


Unter ungeschicktesten Bedingungen würde die Erfassung genau vor dem Zeitpunkt des 
Eintretens neuer Individuen in eine Population erfolgen. Die zuletzt neu eingetretenen 
Individuen würden somit der vollen Mortalität bzw. Emigration unterliegen (Gleichung 5). 


Die Zahl der zwischen i und i+] neu in die Population eingetreten Individuen errechnet 
sich (analog Gleichung 8) für diesen Fall: 
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y 
B; = B, *(()/02)*4z) 07 
j=0 
Und der gesamte Zuwachs an Individuen zwischen 1 und 1+1 schlieBlich ist: 


Big, =B +A y (18) 
l-ġ; 6.7 


l-0 Y. 
1 


Der Term wird Korrekturfaktor Ko; genannt. 


Erfassung und Eintritt von Individuen erfolgen in zeitlichem Abstand 


In den meisten Fallen wird man aber weder direkt nach noch direkt vor dem letzten 
Schlupf der Individuen einen Fang-Wiederfang-Durchgang durchführen, sondern der 
Fang-Wiederfang-Durchgang wird in einem gewissen Abstand t nach dem 
Schlupfzeitpunkt erfolgen (Abb. 1). Nach dem letzten Schlupf vor dem Fang- 
Wiederfang-Durchgang wird also auch der Verlust an Individuen bis zum 
Erfassungszeitpunkt nur den Anteil t an der Überlebenswahrscheinlichkeit #7 betragen, 
die Individuenzahl also um den Faktor ®;', reduziert. 


Anzahl neu In die Populatlon eintretender Individuen 





Tage zwischen i-ter und (H1}-ter Erfassung 


Abb. 1: Die Summe des Zuwachses an Individuen zu den Eintretungsereignissen Z1, Z2, ..,Z5 
addieren sich zu dem Gesamtzuwachs Bices zwischen den zwei Erfassungsereignissen 1 und ı+1 auf. 
Da aber die neu hinzukommenden Individuen einer Mortalität unterliegen, beträgt die zum Erfas- 
sungsereignis it] vorgefundene Zahl an Individuen nur noch Bj. 
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Analog Gleichung 5 ergibt sich: 
B=((...((Bz*®z+Bz)*®z+Bz)*®z+...+Bz)*®,)+Bz)*®z (19) 
Die Gleichung wird aufgelöst: 


B, =B, *(9'6))*45) CO) 


j=0 





ee ee 
Es errechnet sich: Bige = B; 


< 





und KI, = a a wird spezieller Korrekturfaktor genannt. 
1-$; p? 

Im vorliegenden Fall wird davon ausgegangen, dass sich die Verweilwahrscheinlichkeit 
nicht im Tagesverlauf ändert. Dies ist eine vereinfachte Annahme, die meistens streng 
genommen nicht zutreffen wird, denn Emigration findet zu bestimmten Tageszeiten statt 
und Mortalität wird ebenso zu gewissen Zeiten erhöht sein. In den seltensten Fällen wird 
indes ein funktioneller Zusammenhang bekannt sein, wie sich die Verweilwahrschein- 
lichkeit im Tagesverlauf ändert. Daher scheint die Beibehaltung der getroffenen An- 
nahme gerechtfertigt, auch schon dahingehend, um das Berechnungsverfahren nicht über 
Gebühr zu komplizieren. 


Es sei noch darauf hingewiesen, dass bei den zuvor besprochenen optimistischen bzw. 
pessimistischen Erfassungsbedingungen (t = 0 bzw. 1) die Annahme einer im Tages- 
verlauf konstant bleibenden Verweilwahrscheinlichkeit nicht getroffen werden muss. 


Verhalten der Korrekturfaktoren 


Somit sind die Korrekturfaktoren unter pessimistischen (Ko;) und idealen (Ku;) Bedin- 
gungen lediglich Spezialfälle fiir t = 1 bzw. t = 0 und stellen die obere und untere Grenze 
dar, zwischen denen sich Kl;bewegt (Abb. 2). 


Für die Berechnung des Zuwachses zwischen i und i+1 ergeben sich gravierende Unter- 
schiede, von welcher Zeit t zwischen Eintritt von Individuen und Erfassungsdurchgang 
ausgegangen wird. Besonders kleine Verweilwahrscheinlichkeiten ®; wirken sich stark 
auf die Größe es Korrekturfaktors aus, wie Abb. 2 verdeutlicht. Mit zunehmender 
Anzahl von Tagen zwischen zwei Erfassungen nähern sich die obere und untere Grenze 
des Korrekturfaktors an, womit unterschiedliche Zeitspannen t zwischen Eintreten der 
Individuen in die Population am Tag y nicht mehr zu gravierenden Unterschieden führen. 


Wie aber nochmals verdeutlicht wird, führt das unkorrigierte Schätzverfahren nach 
SETTELE et al. (2000) (sozusagen eine Schätzung mit Korrekturfaktor = 1) selbst unter 
der optimistischen Annahme, dass alle am Erfassungstag neu hinzugekommenen Indivi- 
duen in die Stichprobe eingehen können, zu einer Unterschätzung der Gesamtpopula- 
tionsgröße, sobald mehr als ein Tag zwischen i und i+1 liegen und 9; < 1 ist. 
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Tage y zwischen i und i+1 


Abb. 2: Untere Grenze Ku; (jeweils untere Reihe einer Markierung) und obere Grenze Ko; (jeweils 
obere Reihe) der Korrekturfaktors K1; für ausgewählte Überlebenswahrscheinlichkeiten ®i (0,16; 
0,32 und 0,64). Die Geraden geben den Korrekturfaktor K2; an, welcher der Grenzwert von Kl; für 
y> ooist. 


Allgemeiner Korrekturfaktor 


Bei dem oben genannten Korrekturverfahren wurde davon ausgegangen, dass Individuen 
zu bekannten diskreten Zeitpunkten neu in die Population eintreten. 


Was ist aber, wenn sich das Eintreten von Individuen über die gesamt Zeitspanne zwi- 
schen zwei Erfassungsereignissen i und i+1 erstreckt? Um ein möglichst einfaches ma- 
thematisches Modell zu entwickeln, wird angenommen, dass jedes Eintrittsereignis zu 
jeder Zeit zwischen i und i+} mit der gleichen Wahrscheinlichkeit eintreten kann. 
Können zu jeder Zeit Individuen neu in die Population eintreten, so geht die Zeitspanne 
zwischen zwei Eintrittsereignissen j und j+! gegen 0. Das ist das Selbe, wie wenn sich 
die Zahl der Zeitpunkte y, an denen Individuen in die Population eintreten können, ins 
Unendliche erhöht. Angewendet auf Gleichung (15) bedeutet dies: 








_ ay 
Boa = B, "lijm *y (22) 
aT, 
1 
_AY 
Bice = B; * l ia 9; 
1-9, ym 


y 
Es wird die Regel von de l'Hospital (BRONSTEIN et al. 2000) 
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Fa 
timl 2 | = im ) angewendet. Die Differentiale ergeben sich 


folgendermaßen: 


1 


1 1 z “Ay 
Lye) und E Nyon. 
y y y 


1 
! | =1n,*9/ 
2 


í ~ 9 dy y t i 
u lim) 5 lim 2 lim[-1 : 2 = -In6;,]1m®r ° 
yoo ae dy y> — yo» yor 
9 
t 
Da |i (9; ) = 1, bleibt lediglich -Ind, und Gleichung (22) kann vereinfacht werden zu 
y-ız 
-In 6, 
B.... = B, *——+ 23 
iGes i pa db, ( ) 


Wird der Limes aus Gleichung (18) errechnet, führt dies selbstverstandlich zum selben 
Ergebnis. 


Ing; 


Hiermit liegt ein Korrekturfaktor K2, = vor, der unabhängig von der Dauer 


zwischen zwei Erfassungszeitpunkten 1 und i+] ist, sondern lediglich von der Verweil- 
wahrscheinlichkeit zwischen diesen beiden Erfassungszeitpunkten abhängt. 


-In(phi)/(1 - phi) 


K2 
4.5 


0.25 0.5 0.75 
phi 
Abb. 3: Abhängigkeit des Korrekturfaktors K2; von der Verweilwahrscheinlichkeit #i. 
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Wie aus Abb. 3 ersichtlich ist, steigt für kleiner werdende Verweilwahrscheinlichkeiten 
der Korrekturfaktor stark an. Während oberhalb eines ®; > 0,35 der Graph weitgehend 
linear abfällt, macht sich unterhalb von ®; < 0,35 stark der exponentielle Anstieg be- 
merkbar. Geringfügig verschieden geschätzte Verweilwahrscheinlichkeiten wirken sich 
unterhalb eines $; < 0,35 damit deutlich auf den geschätzten Zuwachs der Populations- 
größen B; aus und erhöhen daher stark die Unsicherheit der Schätzung. Verweilwahr- 
scheinlichkeiten $; < 0,35 sollten deshalb vermieden werden. Durch kürzere Abstände 
zwischen zwei Erfassungszeitpunkten i und ı+1 lässt sich ®; größer halten. 


Was leisten die Korrekturfaktoren (nicht)? 


Die oben vorgestellte Methode berechnet den Zuwachs an Individuen zwischen i und 1+] 
bei folgenden Annahmen: Zwischen 1 und 1+] ist der Zuwachs an Individuen und die 
Verweilwahrscheinlichkeit zu jedem Zeitpunkt bzw. zu jedem diskreten Eintretungs- 
ereignis gleich groß. Diese Annahmen werden getroffen, weil zwischen 1 und i+1 keine 
besseren Aussagen über Bz; und $; möglich sind. Tatsächlich variieren aber Bz; und ®,, 
zwischen zwei Erfassungsereignissen. 


B; ist lediglich das Ergebnis aus _ ((..((Bz,*%2,+Bz.)*$7.+Bz;)®25t... 
+Bzy.1)*® zy1)+Bz,)*®;. War Bz; anfangs höher, liefert die Gleichung für B; einen niedri- 
geren Wert als wenn Bz; zu einem späteren Zeitpunkt höher war, selbst wenn ®,; zu allen 


y 
Zeitpunkten konstant war und B..., = > B; in beiden Fällen gleich groß war. ®; schließ- 


jal 


lich ist IIe, . War beispielsweise $7; anfangs höher und später niedriger, so ergibt sich 
j=l 

für B; ein höherer Wert, als wenn ungekehrt 67; anfangs niedriger und später höher war, 

selbst wenn B,, unverändert blieb. In beiden Fällen aber ist ®; gleich groß. Es gibt viele 

Möglichkeiten B,; und #7; zu variieren und dabei zu gleichen B; und ®; zu gelangen. 

Dadurch stellt sich ein bisher noch nicht gelöstes Problem, Konfidenzintervalle für die 

Schätzung der Gesamtpopulationsgröße angeben zu können. 


Die Konfidenzintervalle nach SETTELE et al. (2000) für Gesamtpopulationsgrößen beru- 
hen auf der Annahme, dass B;=B;,., (also neu zwischen i und i+1 hinzukommende Indı- 
viduen keiner Mortalität unterliegen) und berücksichtigen nicht oben genannte Unsicher- 
heiten. Diese Konfidenzintervalle für Gesamtpopulationsgrößen sind daher zu eng. Sie 
können zumindest aber solange eine Orientierungshilfe sein, bis ein besseres Verfahren 
zu Abschätzung der Konfidenzintervalle von Gesamtpopulationsgrößen entwickelt wor- 
den ist. 


Die unkorrigierte Schätzung der Gesamtpopulationsgröße nach SETTELE et al. (2000) 
führt zwangsläufig zu einer Unterschätzung des Zuwachses Bice. Mit Hilfe des speziel- 
len oder des allgemeinen Korrekturfaktors ist die Berechnung eines verbesserten Schätz- 
wertes des Zuwachses und somit der Gesamtpopulationsgröße möglich. Die Korrektur- 
faktoren vermögen also lediglich den dem bisherigen Verfahren anhaftenden Bias, der 
zwangsläufig zu einer Unterschätzung der Populationsgröße führt, zu verkleinem. Die 
Konfidenzintervalle selbst jedoch lassen sich nicht enger fassen. Die vorgestellten 
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Methoden erlauben also eine akkuratere Schätzung, während aber die Präzision nicht 
erhöht wird. 


Empfehlungen zur Durchführung einer Gesamtpopulationsgrößenschätzung 


Die Präzision von Schätzungen von Tagespopulationsgrößen steigert sich mit zunehmen- 
der Intensität (zunehmendem Anteil markierter und wieder gefangener Tiere) einer Fang- 
Wiederfang-Studie. Die Schätzungen von Gesamtpopulationsgrößen verbessern sich mit 
zunehmender Häufigkeit von Erfassungen. Mit anderen Worten ausgedrückt, die Schät- 
zung verbessert sich, indem man die Verweilwahrscheinlichkeit ®; zwischen i und i+1 
erhöht, was durch kürzere Abstände zwischen zwei Erfassungsdurchgängen realisiert 
wird. 


Dem Wunsch nach höherer Präzision steht der damit verbundene höhere Arbeitsaufwand 
entgegen. Daher muss ein Kompromiss zwischen Anzahl und Intensität der Begehungen 
und Präzision geschlossen werden. 


Bei der Planung und Durchführung einer Gesamtpopulationsgrößenschätzung sollten 
folgende Punkte bedacht werden: 


Erwartete Verweilwahrscheinlichkeit ©, 


Wie Abb. 3 zeigt, bedarf es mit sinkender Verweilwahrscheinlichkeit zwischen i und 1+1 
einer ungefähr potentiell ansteigenden Korrektur bei der Berechnung des Zuwachses. 
Unterhalb einer Verweilwahrscheinlichkeit von 0,35 zwischen i und i+] verlässt der 
Graph der Funktion von K2; den annähernd linearen Bereich. Bei kleineren ®; führt eine 
nur geringfügige Unsicherheit in der Schätzung der Verweilwahrscheinlichkeit — wie sie 
immer auftritt - zu einer relativ großen Unsicherheit in der Größe des Korrekturfaktors 
und somit in der Größe des Zuwachses Bige. Soll eine Minimierung des Arbeitsaufwan- 
des angestrebt werden, so sollte trotzdem ®; nicht kleiner als 0,35 angestrebt werden. 


Für Populationen, bei denen eine altersunabhängige Mortalität gleichzeitig anwesender 
Individuen vorliegt, was bei Populationen von Tagfalterarten meistens der Fall sein 
dürfte, gilt folgender Zusammenhang (PFEIFER ct al. 2000) 

6,=$, (24) 

Kenn man die ungefähre Tagesverweilwahrscheinlichkeit ®, der untersuchten Art, so 


lässt sich die Anzahl an Tagen, die bei vorgegebenem #; zwischen zwei Erfassungen 
liegen, errechnet: 


_ Ind; 
eae n (25) 


Wird wert auf eine präzisere Schätzung der Gesamtpopulationsgröße gelegt werden, so 
muss selbstverständlich ein hoher Wert für $; vorgegeben werden. 


Erwartete Tagespopulationsgröße 


Allerdings sollte ein $; von 0,35 nur bei sehr großen Tagespopulationen und eines 
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relativen hohen Anteils markierter Individuen angestrebt werden, da man sonst in Gefahr 
läuft aufgrund von Stochastizität zwischen zwei Erfassungen nur wenige Wiederfänge zu 
erzielen und die Konfidenzintervalle der Jolly-Seber-Schätzung sehr breit werden. 


Ebenfalls ist erst bei sehr großen Populationen gegeben, dass aus theoretischen Uberle- 
gungen angenommen werden kann, dass der relative Zuwachs an Individuen über meh- 
rere Tage als konstant angesehen werden kann (Gesetz der großen Zahl). D. h., der Ein- 
fluss der Stochastizität wirkt sich bei sehr kleinen Populationen sehr viel stärker auf die 
Schwankung der Anzahl täglich neu hinzukommender Individuen aus. 


Bei kleinen Populationen oder wenn der Anteil markierter Tiere gering ist, sollte daher 
eine viel höhere Verweilwahrscheinlichkeit zwischen zwei Begehungen angestrebt wer- 
den, indem die Begehungen ın kürzeren Abständen aufeinander folgen. 


Konstanz des Zuwachses B, zwischen i und i+1 


Die vorgestellte Methode basıert auf der Annahme eines konstanten Zuwachses zwi- 
schen i und 1+1. Die Präzision einer Schätzung hängt wesentlich davon ab, wie konstant 
der tägliche Zuwachs an Individuen erfolgte. In der Regel wird der Zuwachs nicht kon- 
stant über eine längere Zeit erfolgen. Die Schlupfrate an Tagfaltern könnte z. B. wit- 
terungsbeeinflusst schwanken. Oder zu Beginn ıst zunächst ein starker Anstieg des Zu- 
wachses zu erwarten. Durch flexible Wahl der Abstände zwischen i und 1+1 sollten die 
Abschnitte zwischen den Begehungen so gewählt werden, dass die Zuwächse B, zwi- 
schen zwei Erfassungszeitpunkten 1 und 1+1 als annähern konstant angenommen werden 
können. 


Konstanz der Verweilwahrscheinlichkeit ®, 


Unter Umständen können während der Flugzeit drastische Änderungen in der Verweil- 
wahrscheinlichkeit auftreten. Beispielsweise könnte sie sich aufgrund eines Kälteein- 
bruchs oder einer Regenperiode erniedrigen. Sind solche Effekte bekannt oder zu erwar- 
ten, können sie bei der praktischen Durchführung der Gesamtpopulationsgrößenschät- 
zung berücksichtigt werden. Es sollte nach Möglichkeit vermieden werden, dass sich 
innerhalb eines Intervalls $, stark verändern. Dies kann geschehen, indem während der 
Durchführung einer Gesamtpopulationsgrößenschätzung prognostiziert wırd (z. Bsp. mit 
Hilfe von Wetterprognosen), ob sich der Parameter ®,;, während eines geplanten 
Intervalls (vorgesehene Zeit zwischen 1 und i+1) voraussichtlich ändern wird. Ist dies der 
Fall, sollte nach Möglichkeit die Intervalllänge geändert, in der Regel verkürzt werden. 


Begehung möglichst zeitnah nach dem Schlupf 


Treten Individuen zu bestimmten Zeitpunkten in die Population ein, wird der Korrektur- 
faktor K1; angewendet. Findet die Erfassungsbegehung nicht unmittelbar vor oder nach 
dem Schlupf statt, findet bei der Berechnung des Korrekturfaktors der Zeitabstand t 
Eingang in die Formel. Dabei wird von der Annahme ausgegangen, dass die Verweil- 
wahrscheinlichkeit sich innerhalb eines Tages nicht ändert. Diese Annahme wird jedoch 
unter Freilandbedingungen mehr oder weniger stark verletzt. Wenn es daher möglich ist, 
sollte eine Begehung möglichst zeitnah zu dem Eintreten von Individuen in die Popula- 
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tion erfolgen, da dann der Einfluss im Tagesverlauf schwankender Verweilwahrschein- 
lichkeiten gering ist. Die Nähe zum Eintrittszeitpunkt sollte zudem nach diesem gewählt 
werden, da zu diesem Zeitpunkt die Anzahl anwesender, neu in die Population eingetre- 
tener Individuen am Größten sein wird. 


Begehungen am Anfang und Ende der Flugzeit 


Eine Schätzung nach Jolly-Seber ist aus methodischen Gründen erst ab der zweiten Be- 
gehung möglich. (Gleichung 3). Ebenso geht der bei der letzten Begehung erfolgte 
Zuwachs nicht mehr in die Schätzung ein. Wenn die Gesamtpopulationsgröße über eine 
ganze Flugzeit hinweg ermittelt werden soll, dürfte letzteres Problem ın den meisten 
Fällen weniger Bedeutung haben, denn zum Ende eine Flugperiode kommen höchstens 
noch wenige Individuen hinzu. Zum Beginn einer Phänologie kann aber bei den ersten 
Begehungen bereits eine bedeutende Tagespopulationsgröße vorliegen, die zudem einen 
starken Zuwachs erfährt (PFEIFER et al. 2000). Um den Startwert N, zu Beginn einer Ge- 
samtpopulationsgrößenschätzung möglichst frühzeitig berechnen zu können, ist es rat- 
sam, die erste und zweite Begehung möglichst kurz hintereinander durchzuführen. 
Eventuell ist dies durch zwei Begehungen an einem Tag zu verwirklichen. Vorrausset- 
zung ist allerdings, dass die bei der ersten Begehung markierten Tiere sich bis zur zwei- 
ten Begehung gut unter die unmarkierten Tiere gemischt haben. 


Berechnung der Gesamtpopulationsgröße 


Zunächst werden die erforderlichen Parameter zu den Erfassungszeitpunkten mit Hilfe 
eines Software-Programms, welche Berechnungen mit der Jolly-Seber-Methode durch- 
zuführen ermöglicht (zum Beispiel Krebswin (KREBS 1998)), geschätzt. Man erhält die 
Größen N, und B; und %; zu den Erfassungszeitpunkten. Wird davon ausgegangen, dass 
Individuen zu konkreten Zeitpunkten in die Population eintreten, werden mit Hilfe der 
Tage y; zwischen zwei Erfassungszeitpunkten, den Zeiten t; zwischen dem letzten Ein- 
tretungsereignis von Individuen vor einer Erfassung und den Verweilwahrscheinlich- 
keiten ®; die Korrekturfaktoren K1; berechnet. Erfolgt der Eintritt von Individuen mehr 
oder weniger kontinuierlich, werden zur Berechnung der Korrekturfaktoren K2; lediglich 
die Verweilwahrscheinlichkeiten ®; benötigt. Die korrigierten Zuwächse Bices werden 
durch Multiplikation von B; mit K1; bzw. K2; berechnet. Die Gesamtpopulationsgröße 


k-2 
wird mit der Formel Nges = N; + > Bices berechnet. 
i=2 
Ausblick 


Die vorgestellte Methode zur Berechnung einer Gesamtpopulationsgröße ist eine 
Verbesserung der Methode nach SETTELE et al. (2000), indem sie berücksichtigt, dass 
zwischen zwei Erfassungsperioden neu hinzukommende Individuen ebenso wie der Rest 
der Population einer Mortalität unterliegen. Die vorgestellte Methode selbst ist aber 
ebenfalls noch weiter verbesserungsfähig. Absehbare Verbesserungen wären: 


127 


e Eine verbesserte Methode zur Berechnung von Konfidenzintervallen, die berück- 
sichtigt, dass der Zuwachs B; an den einzelnen Eintretungsereignissen zwischen i 
und i+1 nicht konstant erfolgt, sondern stochastischen Prozessen unterliegt. 


e Berücksichtigung, dass der Zuwachs B,; an den einzelnen Eintretungsereignissen 
zwischen 1 und i+1 nicht konstant bleibt, sondern einer Funktion folgt, die im Nach- 
hinein aus der Verlauf der festgestellten Zuwächse zwischen den Erfassungsereig- 
nissen erstellt wird. 


Zusammenfassung 


Bei den meisten Tagfalterarten sind während ihrer Phänologie niemals alle Individuen einer Popu- 
lation gleichzeitig anwesend. Darin unterscheiden sie sich von vielen sonst mit Fang-Wiederfang- 
Methoden analysıerten Tierarten. Die Schätzung einer Gesamtpopulationsgröße bei Tagfaltern 
stellt daher besondere methodische Anforderungen. 


Eine geeignete Schätzmethode für Gesamtpopulationsgrößen wird ausgehend von der Jolly-Seber- 
Methode entwickelt. Mit Hilfe der Jolly-Seber-Methode lässt sich die Populationsgröße N; zum 
Zeitpunkt der Datenerhebung schätzen. Außerdem wird der Individuenzuwachs B; am nächsten 
Erfassungszeitpunkt i+1 und die Verweil(-Überlebens)-Wahrscheinlichkeit ®; zwischen i und i+] 
geschätzt. SETTELE et al. (2000) stellten folgendes Berechnungsverfahren für Gesamtpopulations- 


k-2 
größen auf: N.„=N,+ $B, . Es werden also alle Zuwächse über den gesamten 
i=? 
Erfassungszeitraum aufsummiert. Diese Berechnungsverfahren berücksichtigt aber nicht, dass 
zwischen zwei Erfassungszeitpunkten i und i+} neu hinzukommende Individuen selbst auch einer 
Mortalität unterliegen oder aus der Population emigrieren. Dadurch unterschätzt die Berechnung 
nach SETTELE et al. (2000) die Gesamtpopulationsgröße um die zwischen i und i+1 zwar in die 
Population eingetretenen aber auch wieder verloren gegangenen Individuen. 


In dieser Arbeit wurden daher Korrekturfaktoren für B; entwickelt. 


Für den Fall, dass Individuen zu festen Zeitpunkten in die Population eintreten — etwa bei Tagfal- 


terarten, die zu festgelegten Zeiten Schlüpfen - und dazwischen Intervalle ohne Eintritt von Indivi- 
i 


-0 
duen liegen, kann die Schätzung mit Hilfe des speziellen Korrekturfaktors KI; PEA $, SA 


t 
1-6, 5 
verbessert werden. Wobei y ist die Anzahl an Tagen zwischen 1 und i+] ist und t die Zeitspanne (in 
Bruchteil eines Tages) zwischen Schlupf und Erfassungsdurchgang. 


K2, = un: ist der Korrekturfaktor für Arten, die Jederzeit in die Population eintreten können. Er 


sollte auch in den Fällen herangezogen werden, bei denen der Zeitpunkt des Eintretens neuer 
Individuen nicht bekannt ist. 


Der zeitliche Abstand zwischen Begehungen stellt ein Kompromiss dar zwischen dem Wunsch 
nach hoher Präzision der Schätzung und dem Wunsch den Arbeitsaufwand, d. h. die Anzahl der 
Begehungen während einer Flugperiode, möglichst gering zu halten. Bei einer kleinen Uberle- 
benswahrscheinlichkeit der untersuchten Art darf der Abstand zwischen zwei Begehungen nur kurz 
sein, ebenso wenn die Schätzung präzis sein soll. Außerdem ist die zu erwartende Tagespopula- 
tionsgröße zu berücksichtigen. Ist diese klein, so ist ebenfalls der Abstand zwischen zwei Bege- 
hungen kurz zu halten. 
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